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Abstract. We give an elaborated treatment of discrete isothermic surfaces and 
their analogs in different geometries (projective, Mobius, Laguerre, Lie). We find 
the core of the theory to be a novel projective characterization of discrete isothermic 
nets as Moutard nets. The latter belong to projective geometry and are nets with 
planar faces defined through a five-point property: a vertex and its four diagonal 
neighbors span a three dimensional space. Analytically this property is equivalent to 
the existence of representatives in the space of homogeneous coordinates satisfying 
the discrete Moutard equation. Restricting the projective theory to quadrics, we 
obtain Moutard nets in sphere geometries. 

In particular, Moutard nets in Mobius geometry are shown to coincide with 
discrete isothermic nets. The five-point property in this particular case says that 
a vertex and its four diagonal neighbors lie on a common sphere, which is a novel 
characterization of discrete isothermic surfaces. Discrete Laguerre isothermic sur- 
faces are defined through the corresponding five-plane property which requires that 
a plane and its four diagonal neighbors share a common touching sphere. Equiva- 
lently, Laguerre isothermic surfaces are characterized by having an isothermic Gauss 
map. We conclude with Moutard nets in Lie geometry. 
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1 Introduction 



This paper is a sequel to our paper "On organizing principles of discrete 
differential geometry. Geometry of spheres" IBS21 , where the following dis- 
cretization principles have been formulated: 

• Transformation group principle: smooth geometric objects and their 
discretizations belong to the same geometry, i.e. are invariant with 
respect to the same transformation group. 

• Consistency principle: discretizations of smooth parametrized geome- 
tries can be extended to multidimensional consistent nets. 

Being applied to discretization of curvature line parametrizations of general 
surfaces, these principles led to the definition of principal contact element 
nets. These are nets of contact elements with the property that neighboring 
contact elements share a common sphere. In particular, it was shown that 
the points and the planes of principal contact element nets build circular 
and conical nets, respectively. 

In the present paper we turn to isothermic surfaces which is a special 
class of surfaces admitting a conformal curvature line parametrization. This 
important class of surfaces has been studied by classics |Daj . In the 1990-s 
a relation to the theory of integrable systems has been discovered CGS, 
BHPP| FBPlj . An overview of the modern theory of isothermic surfaces can 
be found in |HJj . The theory has been extended for isothermic surfaces in 
spaces of arbitrary dimension in jBu[ ISchj . In |BPlj the theory has been 
discretized: discrete isothermic surfaces were defined as special circular nets 
with factorized cross-ratios of elementary quadrilaterals. This property is 
manifestly Mobius- invariant. Moreover, it can be consistently imposed on 
three-dimensional nets |HHP|lBP2j . Thus, discrete isothermic surfaces is an 
instance of geometry satisfying both discretization principles. 

In this paper, we give an elaborated treatment of discrete isothermic 
surfaces and their analogs in different geometries (projective, Mobius, La- 
guerre, Lie), applying the discretization principles systematically. We find 
the core of the theory to be a novel projective characterization of discrete 
isothermic nets as Moutard nets. The latter belong to projective geometry 
and are nets with planar faces / : Z 2 — > WF N defined through a five-point 
property: a vertex and its four diagonal neighbors span a three dimensional 
space (thus, in comparison with a generic net with planar faces, the dimen- 
sion drops by one). Analytically this property is equivalent to the existence 
of representatives in the space of homogeneous coordinates y : Z 2 — > IR^" 1 " 1 
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satisfying the Moutard equation 

t\T2V - y || T2y - ny, 

where Tj is the shift in the i-th coordinate direction. The consistency of the 
five-point property follows from the consistency of the Moutard equation. 

Restricting the projective theory to quadrics, we obtain Moutard nets 
in sphere geometries. In particular, Moutard nets in Mobius geometry are 
shown to coincide with discrete isothermic nets. The five-point property in 
this particular case says that a vertex and its four diagonal neighbors lie on 
a common sphere, which is a novel characterization of discrete isothermic 
surfaces. 

In Laguerre geometry discrete surfaces are maps Z 2 — > {planes in R 3 }. 
The planarity of faces in the Laguerre quadric is equivalent to the conical 
property, while the five-plane property requires that a plane and its four 
diagonal neighbors share a common touching sphere. This is a definition 
of discrete Laguerre isothermic surfaces. Equivalently, discrete Laguerre 
isothermic surfaces are characterized by having an isothermic Gauss map. 
The latter class was independently introduced in |WPj . Smooth Laguerre 
isothermic surfaces have been studied in [El IMNl| IMN2j . 

We conclude with Moutard nets in Lie geometry. The latter are special 
Ribaucour sphere congruences Z 2 — > {spheres in R 3 } with the corresponding 
five-sphere property. A particular case is S-isothermic surfaces |BP21 BHS, 

he]. 

2 Discrete Moutard nets 

In considerations of various nets / : Z 2 — > X, we use the following notational 
conventions: for some fixed u G Z 2 , we write / for f(u), and further fa for 
Tif(u) = f(u + a), and for T ~ 1 f( u) = f(u - e*). Also, we freely use 
notations, definitions and results from [BS2 . 

Definition 1 (Discrete Moutard net) A two-dimensional Q-net f : Z 2 —> 
WP N (N > 4) is called a discrete Moutard net, if for any u £ I? the five 
points f and /±i,±2 in a three-dimensional subspace V C RP^, not con- 
taining some (and then any) of the four points f±\, f±2- 

Thus, the defining condition of a discrete Moutard net deals with four el- 
ementary planar quadrilaterals adjacent to one vertex. As a consequence 
of this definition, all nine vertices of such four quadrilaterals of a discrete 
Moutard net lie in a four-dimensional subspace of RP^. 
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Theorem 2 (Discrete Moutard equation) A discrete Moutard net f : 
Z 2 — > MF N possesses a lift to the homogeneous coordinates space y : Z 2 — ► 
M Ar+1 satisfying the discrete Moutard equation: 

2/12 - y = 012(2/2 - 2/1), (1) 

with some a\2 : Z 2 — > R (it is natural to assign the real numbers a\2 to the 
elementary squares of 1?). 

Proof. We start with the observation that for any Q-net / in MF N it is 
always possible (and almost trivial) to find homogeneous coordinates for the 
four vertices of one elementary quadrilateral satisfying the discrete Moutard 
equation on that quadrilateral. Moreover, one can do this for an arbitrary 
choice of homogeneous coordinates for any two neighboring vertices of the 
quadrilateral. Indeed, consider any homogeneous coordinates /, /1, J2, /12 € 
fljAH-l f or vertices of a planar quadrilateral, connected by a linear relation 

/12 = C21/1 + £12/2 + Puf- 

Keep the representatives y = /, y% = f%, say, and set /12 = PyiVyi and 
f2 = oy2 with a = — C21/C12, then y satisfies the discrete Moutard equation 
(^Q) within one elementary quadrilateral. 

Now, for Q-nets with a special property formulated in Definition ^ this 
construction can be extended to the whole net. Start with arbitrary repre- 
sentatives y, yi, and proceed clockwise around the vertex y. We find consec- 
utively: the representatives y_2, 2/1,-2 which assure the Moutard equation 
on the quadrilateral (y, y%, 2/1,-2, V-2)i then the representatives y_i, y-1,-2 
which assure the Moutard equation on the quadrilateral (y, y~\, y-1,-2, y-%), 
and then the representatives 7/2 > 2 which assure the Moutard equation on 
the quadrilateral (y, y_i, y— 1,2 , 2/2) , see Fig.^ In the remaining quadrilateral 
{y, yii 2/12) 2/2)) the representatives 2/, 2/1 > 2/2 are already fixed on the previous 
steps of the construction, so that we can dispose of the representative 2/12 of 
/12 only. Observe that the point with the representative 2/1 — 2/2 belongs to 
the the plane II C MF N of the quadrilateral (/, /1, /12, f2) (obviously) and to 
the three-dimensional space V C WP N through the points /, fi t -2, f-1,-2, 
7-1,2) because of the equation 

2/1-2/2 = (2/1 - y-z) + (y-2 - 2/-1) + (2/-1 - 2/2) 

= 0(2/1-2 - y) + P(y-i,-2 -y) + i(y-i,2 - y)- 

By the hypothesis of the theorem, the point /i ; 2 lies in the latter space V. 
Therefore, the whole line through / and /1 2 lies in the intersection II D V. 
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Figure 1: Constructing a Moutard representative for a projective Q-net with 
a three-dimensional black cross 



Since N > 4, we conclude that in general position liny is the line through / 
and /i2- Thus, the point with the representative y\ — y2 belongs to this line, 
therefore y% — y% is a linear combination of y and y\2- By a suitable choice 
of representative y\% of /12, we can make y\ — 2/2 proportional to y\2 — y- 
Thus, the construction of representatives satisfying the Moutard equation 
closes up around any vertex. This allows to extend the construction to the 
whole lattice 1? . □ 

Definition H is non-applicable in the case when some, and then all of the 
points f±\, f±2 lie in the three-dimensional space V through /, /±i,±2 5 in 
particular, it cannot be used to define discrete Moutard nets in MP 3 . We 
will show that Theorem [2 remains valid if one defines discrete Moutard nets 
in MP 3 as follows. 

Definition 3 (Discrete Moutard net in MP 3 ) A two-dimensional Q-net 
f : 1? — > MP 3 is called a discrete Moutard net, if for any u € 1? the following 
condition is satisfied: the three planes 

n( up ) = (/,/i 2 ,/_i, 2 ), n( down ) = (/,/ 1 ,_ 2 ,/-i,- 2 ), nW = (/,/!,/_!) 

have a common line £^ . 

Remark 1. It is not difficult to see that in the context of Definition 
n with N > A the requirement of Definition |31 is automatically satisfied. 
Indeed, in this case all nine points /, f±\, f±2 and /±i,±2 he in a four- 
dimensional subspace of PJP^. In this subspace one can consider, along 
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with the three-dimensional subspace V, the three-dimensional subspaces 
y(up) containing the two quadrilaterals (/, fi, f 12 , / 2 ), (/, f-i, f-1,2, h), and 
y (down) containing the quadrilaterals (/, /1, /i_ 2 , /- 2 ), (/, f-l, f-1,-2, f-2). 
Obviously, one has: 

jj(up) _ y(up) p y jj(down) _ y(down) p y jj(1) _ y(up) p y(down)_ 

Generically, three three-dimensional subspaces V, y( up ) and y( down ) of a 
four-dimensional space intersect along a line £W, 

Remark 2. There is an asymmetry between the coordinate directions 1 
and 2 in Definition |31 However, this asymmetry is apparent: the condition 
in Definition 01 is equivalent to the requirement that the three planes 

n (lcft) = (/, /_i,2, f-1,-2), n^ ht ) = (/, / 1)2 , / lr _ 2 ), n( 2 ) = (/, / 2 , /_ 2 ) 

have a common line^ 2 ). One way to see this is to consider a central projec- 
tion of the whole picture from the point / to some plane not containing /. 
In this projection, the planarity of elementary quadrilaterals (f,fi,fij,fj) 
turns into collinearity of the triples of points fi, fj and fy. The traces 
of the planes n^ up \ n^ down ^ and II^ 1 ) on the projection plane are the lines 
(f 12, f-1,2), (fi -2, f-1,-2), and respectively, and the requirement 

of Definition El turns into the requirement for these three lines to meet in 
a point. Similarly, the traces of the planes n( left ), n( nght ) and II^ 2 ) on the 
projection plane are the lines (f -1,2, f-1,-2), (fi,2, fl-2), and (/ 2 ,/_ 2 ), re- 
spectively. The requirement for the latter three lines to meet in a point is 
equivalent to the previous one - this is the statement of the famous Desar- 
gues theorem, see Fig. |21 

Another way to demonstrate the actual symmetry between the coordi- 
nate directions 1 and 2 in Definition 01 is to show that Theorem [2 still holds 
in MP 3 . Indeed, the discrete Moutard equation Q is manifestly symmetric 
with respect to the flip 1 «-> 2. 

Proof of Theorem |51 for N = 3. We start the proof exactly as in the 
general case N > 4. The only thing to be changed is the demonstration 
of the fact that the point with the representative y 2 — y\ lies on the line 
through / and /i 2 . To do this in the present situation, we first observe that, 
due to 

vi - y-i = (yi - y-2) + (y-2 - y-i) = a (2/1 -2 - y) + &{y-i-2 - y), 

the point with the representative y\ — y-i lies in the plane n( down ). Obvi- 
ously, it lies also in II^ 1 ). Therefore, it lies on the line As a consequence 
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Figure 2: Desargues theorem 



of the property of Definition |31 it belongs also to the plane n( up ). Now from 

2/2 — 3/i = (2/2 - y-i) - (yi - y-i) = i{y-i,2 -y)- (2/1 - 2/-1) 

we find that the point with the representative yi — y\ belongs to n( up ), as 
well. Since the point with the representative yi — y\ also belongs (obviously) 
to the plane of the quadrilateral (/, /i, /12, /2), we conclude that it lies in 
the intersection of the latter plane with n( up ) = (/, /12, /-1 2)) which is, in 
the generic case, the line through / and fyi- □ 

Definitions Q |3] are essentially dealing with two-dimensional Q-nets. 
However, the characterization of discrete Moutard nets given in Theorem 
121 opens a way to define multi-dimensional Moutard nets, and, in particular, 
to define transformations of Moutard nets with remarkable permutability 
properties. Namely, it turns out that eq. (^Q) can be posed on multidimen- 
sional lattices. 

Definition 4 (T-net) A map y : 7L m — > R N is called an m- dimensional 
T-net (trapezoidal net), if for any u G Z m and for any pair of indices i ^ j 
there holds the discrete Moutard equation 

Vij ~ V = a>ij(yj - yi), (2) 

with some aij : Z m — > R, in other words, if all elementary quadrilaterals 
(y,yi,yij,yj) are planar and have parallel diagonals. 
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Of course, coefficients a^- have to be skew-symmetric, etjj = —aji. We 
show that three-dimensional T-nets are described by a well-defined three- 
dimensional system. 

Theorem 5 (Elementary hexahedron of a T-net) Given seven points 
y, yi, and yij (1 < i ^ j < 3) in R , such that eq. is satisfied on 
the three quadrilaterals (y,yi,Vij,yj) adjacent to the vertex y, there exists a 
unique point 2/123 such that eq. is satisfied on the three quadrilaterals 
(yi,yij,yi23,Vik) adjacent to the vertex y^. 

Proof. Three equations for the faces of an elementary cube of Z 3 adja- 
cent to 2/123, give: 

nyjk = (1 + (Tia jk )(aij + a ki ))yi - lT,<i l i,.}<i; j y l - (Tia jk )a ki y k . 

They lead to consistent results for 7/123 f° r arbitrary initial data, if and only 
if the following conditions are satisfied: 

1 + (na23)(ai2 + «3i) = -(T2«3i)ai2 = -(r 3 ai2)a 3 i, 
1 + (T2fl3i)(a23 + 012) = -(T 3 ai 2 )a23 = -(Tia 23 )ai2, 
1 + (73ai 2 )(a23 + a 3i) = -(na23)a3i = -(T2a3i)a23- 

These conditions constitute a system of 6 (linear) equations for 3 unknown 
variables T{aj k in terms of the known ones aj k . A direct computation shows 
that this system is not overdetermined but admits a unique solution: 

Tia23 _ T2Q31 _ 73012 1 ^ 

023 031 Ol2 012CJ23 + 023^31 + 031^12 

With TiOj k so defined, eqs. © are fulfilled on all three quadrilaterals adja- 
cent to yi23- □ 

Eqs. (JHJ) represent a well-defined birational map {ajk} i— > {Tiajfc}, which 
can be considered as the fundamental 3D system related to T-nets. It is 
sometimes called the "star-triangle map". 

Theorem El means that the defining condition of T-nets (parallel diago- 
nals of elementary planar quadrilaterals) yields a discrete 3D system with 
fields on vertices taking values in an affine space R^. This system can be 
considered as an admissible reduction of the 3D system describing Q-nets 
in M. N . Indeed, if one has an elementary hexahedron of an affine Q-net 
y : I? — > R N such that its elementary quadrilaterals (y,Vi,yij,yj) have par- 
allel diagonals, then the elementary quadrilaterals (yi,yij,yi23,Vik) have this 
property, as well. To see this, observe that the point 1/123 from Theorem El 
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satisfies the planarity condition, and therefore it has to coincide with the 
unique point defined by planarity of the quadrilaterals (yi,yij,yi23,yik)- 

The 4D consistency of T-nets is a consequence of the analogous property 
of Q-nets, since T-constraint propagates in the construction of a Q-net from 
its coordinate surfaces. On the level of formulas we have for T-nets with 
m > 4 the system (J2J), while the map {a,jk\ i— > {Tia^} is given by 

I^B = I . ( 4 ) 

All indices k vary now between 1 and m, and for any triple of pair- 
wise different indices (i,j,k), equations involving these indices solely, form 
a closed subset. 

The multidimensional consistency of T-nets yields in a usual fashion Dar- 
boux transformations with permutability properties (which in the present 
context should be called discrete Moutard transformations). We refer to 
BS1 for the background on the relation of multidimensional consistency to 
Darboux transformations, and give here only the formulas for the discrete 
Moutard transformation of eq. (j2J into 

vti - y + = a ti(yt - vt)- ( 5 ) 



These formulas read: 



yf ~y = bi{y + -yi), (6) 



where the quantities 6, and the transformed coefficients afj are defined by 
equations 

nty _ a tj _ 1 / 7 n 

bj a {j {hi - bj)uij + bib, ' 

It is not difficult to recognize in eqs. (JSJ), © the same Moutard equations 
PJ) on the (m + l)-dimensional lattice, with the superscript "+" used to 
denote the shift T m +\- Similarly, eqs. (J7J) are nothing but the star-triangle 
formulas ® with 6j = aj i?n+ i. 

Remark. We learned about the projective five-point characterization 
of two-dimensional Moutard nets from conversations with A. Doliwa. It 
should be noted that three- and higher-dimensional Moutard nets admit a 
different projective characterization (planarity of tetrahedra formed by odd 
or by even vertices of any elementary cube), see |Do3j . Also the paper |Do2j 
by A. Doliwa deals with a closely related notion of discrete Koenigs nets. 
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2.1 Discrete Moutard nets in quadrics 

We have seen that discrete Moutard nets (or, more precisely, their T-net 
representatives) constitute an admissible reduction of Q-nets. The restric- 
tion to a quadric constitutes another admissible reduction |Dolj . Imposing 
two admissible reductions simultaneously, one comes to T-nets in quadrics. 
Let M. N be equipped with a non-degenerate symmetric bilinear form (•, •) 
(which does not need to be positive-definite), and let 

oN 



Q = {y Gl' v : (y,y) = k } 



be a quadric in M. . We study T-nets y : Z m — > Q. This leads to a dis- 
crete 2D system, since constructing elementary quadrilaterals of T-nets in Q 
corresponding to elementary squares of the lattice Z m admits a well-posed 
initial value problem: given three points y, yi,?/2 £ Qj one finds a unique 
fourth point yi2 € Q, yi2 7^ 2/, satisfying the discrete Moutard equation 

2/12 - 2/ = 012(2/2 - yi)- 

Indeed, the condition 

(2/12, 2/12) = (y + 012(2/2 - yi),y + ^12(2/2 - yi)} = ^ 

leads to a quadratic equation for a±2, which has one trivial solution a 12 = 
<^ yi2 = V, and one non-trivial: 



a\2 



(y, yi - yi) 

KQ - (1/1,2/2) 



This elementary construction step, i.e., finding the fourth vertex of an ele- 
mentary quadrilateral out of the known three vertices, is symbolically rep- 
resented on Fig. |21 




Figure 3: 2D system on an elementary quadrilateral 
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Turning to an elementary cube of dimension m > 3, we see that one 
can prescribe all points y and y, for all 1 < i < m. Indeed, these data are 
independent, and one can construct all other vertices of an elementary cube 
from these data, provided one does not encounter contradictions. To see the 
possible source of contradictions, consider in detail the case of m = 3. From 
y and (1 < i < 3) one determines all by 



Vij - y = Oij(yj - y%), 



(y,yi - yj) 



(8) 



After that one has, in principle, three different ways to determine yi23, from 
three squares adjacent to this point; see Fig. These three values for yi23 
have to coincide, independently of initial conditions. 

Definition 6 (3D consistency) A 2D system is called 3D consistent, if it 
can be imposed on all two-dimensional faces of an elementary cube ofl?. 




1.23 



Figure 4: 3D consistency of 2D systems 



There holds a quite general theorem, analogous to Theorem 5 of BS2 : 

Theorem 7 (3D consistency yields consistency in all higher dimen- 
sions) Any 3D consistent discrete 2D system is also m-dimensionally con- 
sistent for all m > 3. 

Proof goes by induction in m and is analogous to the proof of Theorem 5 
from |BS2] . □ 

Theorem 8 (T-nets in quadrics are 3D consistent) The 2D system 
10) governing T-nets in Q is 3D consistent. 
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Proof. This can be checked by a tiresome computation, which can be 
however avoided by the following conceptual argument. T-nets in Q are 
a result of imposing two admissible reductions on Q-nets in M N , namely 
the T-reduction and the restriction to a quadric Q. This reduces the effec- 
tive dimension of the system by 1 (allows to determine the fourth vertex 
of an elementary quadrilateral from the three known ones), and transfers 
the original 3D equation into the 3D consistency of the reduced 2D equa- 
tion. Indeed, after finding yi2, 2/23 and 2/13, one can construct 2/123 according 
to the planarity condition (as intersection of three planes). Then both the 
T-condition and the Q-condition are fulfilled for all three quadrilaterals ad- 
jacent to 2/123, according to Theorem |S1 and the result of |Dolj . Therefore, 
these quadrilaterals satisfy our 2D system. □ 

We mention also an important property of T-nets in quadrics used in 
the sequel: the functions 

ai = {y,Vi), (9) 
defined on edges of Z m parallel to the i-th coordinate axes, satisfy 

naj = aj, i^j, (10) 

i.e, any two opposite edges of any elementary square carry the same value 
of the corresponding 04. Indeed, equations 

(yij,yj) = (yi,y), (yij,yi) = (vj,y)- 

follow from (JHJ) by a direct computation. 

3 Isothermic surfaces in Mobius geometry 

Definition 9 (Discrete isothermic surface) A two-dimensional circular 
net f : 1? — > M. N is called a discrete isothermic surface, if the corresponding 
f : 7? — > L, N+1,1 is a discrete Moutard net. 

Prom Definitions ^ 01 there follows a geometric characterization of dis- 
crete isothermic nets: 

Theorem 10 (Central spheres for discrete isothermic nets) 

(a) A circular net f : 1? — > M. N not lying on a two-dimensional sphere 
is a discrete isothermic net, if and only if for any u € Z 2 the five points f 
and /±i,±2 He on a two-dimensional sphere not containing some (and then 
any) of the four points f±\, f±2- 
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(b) A circular net f : 1? — > 8 2 C is a discrete isothermic net, if and 
only if for any u £ Z 2 t/ie three circles through f , 

C (up) = drde(/j /l2) C (down) = drcle(/) /ij 2) / _ i _ a)) 

C« = circle^, /!,/_!), 

/iai>e one additional point in common, which is also equivalent for the three 
circles through f, 

C (left) = drde(/) / ij 2)) fright) = drde(/) /l 2) /i5 _ 2)) 

C( 2 ) = circle(/,/ 2 ,/_ 2 ), 
to have one additional point in common. 



Figure 5: Four circles of a generic discrete isothermic surface, with a central 
sphere. 

The cases a), b) of Theorem ITU1 are illustrated on Figs. ISJIH1 respectively. 

Another characterization of discrete isothermic surfaces can be given 
in terms of the cross-ratios. Recall that for any four concircular points 
/) fi > f2 1 fi2 £ t w their (real- valued) cross-ratio can be defined as 

<?(/, fl, fl2, h) = (fl ~ /)(/l2 " fl)-\fl2 ~ / 2 )(/2 " f)' 1 . (11) 

Here multiplication is interpreted as the Clifford multiplication in the Clif- 
ford algebra Q£(M. N ). Recall that for x, y £ M N the Clifford product satisfies 
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Figure 6: Four circles of a planar discrete isothermic net. 



xy + yx = —2(x, y), and that the inverse element of x € M in the Clifford 
algebra is given by x _1 = —x/\x\ 2 . Alternatively, one can identify the plane 
of the quadrilateral (/, fx, /12, /b) with the complex plane C, and then inter- 
pret multiplication in eq. (jl 1|) as the complex multiplication. An important 
property of the cross-ratio is its invariance under Mdbius transformations. 

Theorem 11 (Four cross-ratios of a discrete isothermic net) A cir- 
cular net f : 1? — > is a discrete isothermic net, if and only if the 
cross-ratios q = q(f , fi, f 12, f 2) of its elementary quadrilaterals satisfy the 
following condition: 

Q " Q-1,-2 = q~i ■ Q-2- (12) 

Here, like in Sect. the negative indices —i are used to denote the backward 
shifts t^ 1 , so that, e.g., q-i = q(f-i, /, /2, f-1,2), see Fig.\]\ 

Proof. Perform a Mobius transformation sending / to 00. Under such a 
transformation, the four adjacent circles through / turn into four straight 
lines f±if±2, containing the corresponding points /±i,±2- Formula (|12j) turns 
into the following relation for the quotients of (directed) lengths: 

I/2/12 1 _ l/l/l, -2| _ |/-2/-l,-2| _ |/-l/-l,2| _ ^ 
|/l2/l I |/l,-2/-2| |/-l,-2/-l| |/-1,2/2| 

If the affine space through the points f±i, f±2 is three-dimesnional, then eq. 
(jlHjl is equivalent to the fact that the four points /±i ±2 he in a plane, which 
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Figure 7: Four adjacent quadrilaterals of a discrete isothermic net: the 
cross-ratios satisfy q ■ q-1-2 = q-i • 9-2, the five points / and /±i,±2 he on 
a sphere 



is a sphere through / = 00. This is a three-dimensional generalization of 
the Menelaus theorem; since this generalization is not very well known, we 
give it, with a proof, in appendix to this section. If, on the contrary, the 
four points f±i, f±2 are co-planar, then we are in the situation of Fig. |2l 
described by the Desargues theorem. Here, we apply the Menelaus theorem 
twice, to the triangle /- 1/2/1 intersected by the line /-i,2/i2> and to the 
triangle /-1/-2/1 intersected by the line f~i-2fi,-2- both lines meet the 
line /-1/1 at the same point if and only if 

I/2/12 1 |/-i/-i, 2 j _ |/-il (1) | = |/- 2 /i,- 2 | |/-i/-i,- 2 j 

|/l2/l|' |/-1,2/ 2 | ' K«/l| |/l,-2/l| ' |/-l,-2/-2|' 

This yields JT3J). □ 

The claim of Theorem \H\ can be re-formulated as follows. 

Corollary 12 (Factorized cross-ratios for a discrete isothermic net) 

A circular net / : 1? — > M. N is a discrete isothermic net, if and only if the 
cross-ratios q = q{f , fi, f 12, f 2) of its elementary quadrilaterals satisfy the 
following condition: 

9(/,/i,/i2,/ 2 ) = -, (14) 

«2 

with some edge functions «j satisfying the labelling property il(J\) . 
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Clearly, functions Oj are denned up to a common constant factor. Actually, 
it was this characterization of discrete isothermic nets that was used as a 
definition in the pioneering paper |BPlj . 

Actually, edge functions on in Corollary ^] admit a nice geometric ex- 
pression. According to Theorem |21 a discrete isothermic net / : Z 2 — > \W N 
can be characterized by the existence of representatives y : 1? — > L, N+1,1 in 
the light cone satisfying the discrete Moutard equation l|T]). We will use for 
these T-net representatives the notation 



y 



"7 = a~ 1 (/ + eo + |/| 2 e o). (15) 



Thus, s _1 denotes the eo-component of the T-net representative y of an 
isothermic net /. The function s : 1? — > R plays the role of a discrete 
metric of the isothermic net /. Now, eq. Q defines functions on edges, 

a l = -2{y,y l )= l -^l, (16) 

SSi 

possessing property ()1U|) : any two opposite edges of any elementary square 
carry the same value of the corresponding a%. 

Theorem 13 (Cross-ratios through discrete metric) Edge functions 
on participating in the factorization \14\ ) of the cross-ratios of elementary 
quadrilaterals of a discrete isothermic net can be defined by eq. Mb)) . 

Proof. Comparing the eo-components in the Moutard equation y\2 — y = 

^12(2/2 — yi), we find: a\2 = {$12 ~ s ~ 1 )/( s 2~ — s i )■ Therefore, we can 
re-write the Moutard equation as 

J_ _ ]_\ ffv2 _t\ = (— - 

S 2 SiJ\Si2 S/ \Si2 s)\S2 S\ 

which is equivalent to 

fi - f _|_ fi2 - fi _ h~ f + /12 - h 

SS\ S\Sl2 SS2 S2S12 

The R^-part of the latter equation, i.e., 

fi — f . fi2 — fx f2 — f . fn — h n n 

1 — 1 , {i-O) 

SSi SlSl2 SS2 S2S12 

can be rewritten with the help of eq. (j!6ft as 

fi~f , fi2 - fi h~ f fi2 - h / in x 

Ji-/r 712 -/if 72 - 7 712 - 72 2 
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In terms of the inversion in the Clifford algebra e£(R N ), this can be presented 
as 

- /r 1 + a 2 (/i2 - fiY 1 = « 2 (/ 2 - JT 1 + «i(/i2 - /a) -1 - (20) 

This latter equation is, in the generic case /12 + / 7^ /1 + /2, equivalent to 
eq. I|14p. It is not quite straightforward to show this equivalence in case of 
non-commutative variables / G (^(R^). But one can identify the plane of 
the quadrilateral (/, /1, /i2> / 2 ) with C, and then eq. (fT§)) is the (complex 
conjugate of) eq. 1)20(1 . where now all variables are commutative (complex 
numbers), and in this case the equivalence to eq. (|14|) is immediate. □ 

T-nets in the light cone L Ar+1,1 are 3D-consistent. This yields also the 
3D-consistency of the cross-ratio equation (|14|) with prescribed labelling on 
of the edges, i.e., of the 2D equation 

q(f,fi,fii,fj) = —- (21) 

Otj 

Both constructions provide us with a well-defined notion of multidimensional 
discrete isothermic nets, and therefore with Darboux transformations of 
discrete isothermic nets with the usual permutability properties. 

We finish this section with the notion of duality for discrete isothermic 
nets. 

Theorem 14 (Dual discrete isothermic net) Let f : Z m — ► R w be a 

discrete isothermic net, with the T-net representatives in the light cone 

y = s- 1 / = s~\f + e + |/| 2 eoo ) : Z m - h N+1,1 . 
Then the W N -valued discrete one- form 8f* defined by 

Sif* = — = a irM\2 ' i = l,...,m, (22) 



is closed. Its integration defines (up to translation) a net f* : Z m — > 
called dual to the net f . The net f* is a discrete isothermic net, with 



q(f*J*,f* j j;) = -- (23) 
a j 

Define also the function s* : Z m — > R as s* = s _1 . Then the net 

y* = (s*)- 1 /* = (s*)- 1 ^* + e + I/* | 2 eoo) : ^ -» L Ar + 1 ' 1 
is a T-net in the light cone. 
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Proof. Clearly, for any pair of indices i,j the function / satisfies an equation 
analogous to eq. Q17|). which expresses the closeness of the M Ar+1 ' 1 -valued 
one- form defined by 6if /(ssj). Unfortunately, the net obtained by integra- 
tion of this one-form does not lie, in general, in the light cone L Ar+1,1 and 
cannot be taken as the dual net /*. We use the following trick for the con- 
struction of the dual net /* in the light cone. The IR^-part of eq. (|17j). i.e., 
eq. pHjl. expresses the closeness of the M^-valued one-form 6{f* = Sif /(ss«), 
whose integration gives a dual net /* in R-^. Eq. (|23j) follows immediately 
from eq. (j22j) and implies that /* is a discrete isothermic net. In particular, 
it is a circular net, so that /* = /* + eo + |/*| 2 eoo is a conjugate net in 
the light cone. It remains to show that the so defined /* is a Moutard net, 
with the T-net representatives y* = (s*) -1 /*. This claim is equivalent to 
closeness of the discrete M Ar+1 ' 1 -valued one-form 5{f* /(s*s*). Since /* is a 
conjugate net in the light cone, it is enough to prove the closeness of the 
IR^-valued one-form 5{f* /(s*s*). But for s* = s^ 1 we have: 

= (ssi)5if* = (ssi)— = 5if, 

S*S* SSi 

which is automatically closed. □ 



3.1 Appendix: generalized Menelaus theorem 

Theorem 15 Let Pq, P n be n + 1 points in general position in W 1 , so 
that the affine space through the points Pi is n- dimensional. Let P^i+x be 
some points on the lines PjPj+i (indices are read modulo n + 1). The n + 1 
points Pa+i lie in an (n — 1) -dimensional affine subspace, if and only if the 
following relation for the quotients of the directed lengths holds: 

TT \PjPj,i+l I _ r_-i\n+X 
l L \Pi,i + lPi + l\ 

Proof. The points Pi^+i lie in an (n — l)-dimensional affine subspace, if 
there is a non-trivial linear dependence 

n n 

J, aiiPij+i = with Oi = 0. 

i=0 i=0 

Substituting Pi.i+i = (l—£i)Pi + £iPi+i, and taking into account the general 
position condition, which can be read as linear independence of the vectors 
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Pi — Pq , we come to a homogeneous system of n + 1 linear equations for n + 1 
coefficients af. 

ZiOti + (1 - = 0, i = 0, . . . , n 

(where indices are understood modulo n + 1). Clearly it admits a non-trivial 
solution if and only if 

4 Isothermic surfaces in Laguerre geometry 

Definitions and constructions of this section can be generalized for higher- 
dimensional nets 

P : Z m -> {hyperplanes in R N }, 

with the most natural geometric case m = N— 1 (discrete Laguerre geometry 
of hypersurfaces in M. N ); we restrict ourselves to the case of surfaces in M 3 , 
i.e., N = 3, m = 2. 

Definition 16 (Discrete L-isothermic surface) A two-dimensional con- 
ical net P : Z 2 — > {planes in M 3 } is called a discrete L-isothermic surface, 
if the corresponding p : H? — > L 4,2 is a discrete Moutard net. 

Recall that for an (oriented) plane P = {x £ M 3 : (v, x) = d} with the unit 
normal vector u£§ 2 and d G R its representative p in the Lie quadric L 4,2 
is given by 

p = v + • e + 2de oa + 1 • e 6 . 

Recall also that the vectors v : Z 2 — > § 2 comprise the Gauss map for a given 
conical net P : Z 2 — > {planes in M 3 }. 

From Definitions Q 01 there follows a geometric characterization of dis- 
crete L-isothermic nets: 

Theorem 17 (Central spheres for discrete L-isothermic nets) 

(a) A conical net P : 1? — > {planes in R 3 } not tangent to a two- 
dimensional sphere is a discrete L-isothermic net, if and only if for any 
u G Z 2 the five planes P and P±i,±2 are tangent to a two-dimensional sphere 
not touching some (and then any) of the four planes P±\, P±2- 
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(b) A conical net P : Z 2 — > {tangent planes of § 2 C R 3 } is a discrete 
L-isothermic net, if and only if for any u £ Z 2 the three cones through P, 

C ( UP ) = cone ( P) P 12 , P_ 12 ), C( down ) = cone(P, P x _ 2 , P_! _ 2 ), 

=cone(P,Pi,P_i), 

Ziaue one additional plane in common, which is also equivalent for the three 
cones through P, 

C (left) = cone(jP) p_ 12) p_ 1 _ 2 ) ) C-(right) = cone(jP) p i 2) p 1 _ 2 ) j 

=cone(P,P 2 ,P_ 2 ), 
to /ia?;e one additional plane in common. 

The (generic) case a) of Theorem El is illustrated on Fig. |HJ 




Figure 8: Five diagonally neighboring planes of a generic discrete L- 
isothermic surface, with a central sphere. 

Theorem 18 (Gauss map of an L-isothermic net is an isothermic 
net in the sphere) Gauss map of an L-isothermic net is a discrete isother- 
mic net in S 2 . Conversely, if for any u € 1? the four planes P, P\, P2, P12 
of a net P : 1? — > {planes in M 3 } meet at a point, and the Gauss map of the 
net P is isothermic, then P is an L-isothermic conical net. 
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Proof. First, let P be an L-isothermic net. Then for some c : Z 2 — > R the 
net cp is a T-net in the Lie quadric. As a consequence, c(y + 1 • eo + 1 • e^) 
is a T-net in the light cone L 4,1 of the Minkowski space M 4,1 of the Mobius 
geometry for N = 3. Therefore, the net v : Z 2 — > 8 2 C R 3 is isothermic. 

Conversely, let the net v : Z 2 — > S 2 C R 3 be isothermic. This is equiva- 
lent to the existence of the function c : Z 2 — > R such that c(v, 1) is a T-net. 
If now (u, x) = d is the equation of the plane P, then the existence of the 
common intersection point of the planes P, Pi, P2, Pi 2 yields that the func- 
tion cc/ satisfies the same Moutard equation as the function cv. Therefore, 
c(v, d, 1) is a T-net, so that p is a discrete Moutard net. □ 

5 Lie geometry: S-isothermic nets 

Two-dimensional nets in the Lie quadric L 4,2 are discrete congruences of 
spheres. An interesting class of such congruences is constituted by discrete 
Moutard nets in P(L 4 ' 2 ). We leave a general study of this class for a future 
research, and describe here as an example a particularly interesting subclass, 
for which the T-net representatives in L 4,2 have a fixed eg-component: 



Omitting the constant and therefore non-interesting eg-component, we come 
to a T-net in a hyperboloid of the Lorentz space of the Mobius geometry, 





{f€R 4 - 1 :<f,0=« 2 }- 



Definition 19 (S-isothermic net) A map 



S : Z 2 — > {oriented spheres in R 3 



} 



is called an S-isothermic net, if the corresponding map 



s : Z- 





(24) 



is a T-net. 



Thus, S-isothermic nets are governed by equation 



I12 - s = ai 2 (s 2 - Si) 



(Mi - S2} 



OL\ - Ct2 



(25) 



an = 



,2 -(Si,S 2 ) k : 



2 - (Si,S 2 ) ' 



21 



with the quantities on = (s, §i) depending on m only. If (oriented) radii of 
all hyperspheres become uniformly small, r(u) ~ ks(u), k — > 0, then in the 
limit we recover discrete isothermic nets. 

Consistency of T-nets in L 4 / (which is a particular case of Theorem 
|SJ) yields, in particular, Darboux transformations for S-isothermic nets. A 
Darboux transform s + : Z m — > L 4 / of a given S-isothermic net s : Z m — > L 4 / 
is uniquely specified by a choice of one of its spheres s + (0). 

We turn now to geometric properties of S-isothermic nets. First of all, 
S-isothermic nets form a subclass of discrete R-congruences of spheres (see 
BS2 for a geometric characterization of discrete R-congruences). Further, 
consider the quantities (s, Sj) which have the meaning of cosines of the inter- 
section angles of the neighboring spheres (resp., of their so called inversive 
distances if they do not intersect). Then these quantities (s,s?) have the 
labelling propery, i.e., depend on m only. 

There holds the following generalization of Theorem 1141 

Theorem 20 (Dual S-isothermic net) Let 

S : Z m -» {oriented spheres in R 3 } 



be an S-isothermic net. Denote the 
of S by c : Z m — > R 3 and r : Z m - 
discrete one-form 5c* defined by 

die = — , 

rri 



Euclidean centers and (oriented) radii 
-» R, respectively. Then the R 3 -valued 

1 < i < m, (26) 



is closed, so that its integration defines (up to a translation) a function 
c* : Z m -> R 3 . Define also r* : 7L m -» R by r* = r _1 . Then the spheres S* 
with the centers c* and radii r* form an S-isothermic net, called dual to S. 

Proof. Consider equation 

&ij <5 — CLij^Sj 5^), (^^) 

in terms of s from (|24f) . Its eo-part yields: = (r^ 1 — r~ l )/(rj l — r^ 1 ). 
This allows us to rewrite eq. ()27(l as 

( r -l _ r - l ){s i3 -§) = (r^ 1 - r- 1 )^ - s t ). (28) 

A direct computation shows that the R 3 -part of this equation can be rewrit- 
ten as 

C i ~ C , C ij ~ C » _ C j ~ C , C U ~ C j z-29) 

' ' I 1 % 1 I J ' ' J 1 J ' % J 
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which is equivalent to closeness of the form 5c* defined by (|26j) . In the same 
way, the e^-part of eq. Q28JI is equivalent to closeness of the discrete form 
8w defined by 

^(|c| 2 -r 2 ) 
biW = , 1 < i < m. 

For similar reasons, the second claim of the theorem is equivalent to closeness 
of the form 

8i(\c*\ 2 - (r*) 2 ) 

biW = — , 1 < % < m, 

r 

where, recall, r* = 1/r. With the help of c* — c* = (q — c)/rrj, one easily 
checks that the forms 5w and 5w* can be written as 

5iW = (c* - c*,Ci + c) - — + — , 

r ri 

5iW* = (at - c,c* + c*) h — . 

rj r 

The sum of these one- forms is closed: 

5 i (w + w*) = 2(c*,c i ) -2(c*,c), 

therefore they are closed simultaneously. □ 

An interesting particular case of S-isothermic surfaces is characterized 
by touching of any pair of neighboring spheres. In this case the limit of 
small spheres is not relevant, therefore it is convenient to restrict the con- 
siderations to a fixed value of k = 1. Clearly, in this case both oti = (s,Ti§), 
i = 1,2, can, in principle, take values ±1. However, it is easily seen from 
()25|) that in case ol\ = ct2 one gets only trivial nets. Thus, we assume that 

{§, §i) = (h, S12) = -1, {s, h) = (h, S12) = 1- (30) 

Interestingly, these touching conditions already yield that the linear depen- 
dence of the spheres has the Moutard shape. 

Theorem 21 S-isothermic surfaces with touching spheres can be character- 
ized by any of the following equivalent descriptions: 

• Q-congruence of spheres ( Q-net in the Lorentz space of Mobius geom- 
etry) with touching spheres 

• R-congruence (Q-net in Lie quadric) with touching spheres 

• T-net in the Lie quadric with touching spheres 



23 



Proof. Let s, s\, s 2; §12 be four oriented spheres in L 1 ' , pairwise touching 
so that eq. (|3U|) is fulfilled, and linearly dependent. (We remark that in the 
present situation the geometric meaning of linear dependence is the existence 
of a common orthogonal circle through the touching points.) We show that 
the linear dependence has to be of the Moutard form: 

S12 - s = Oi 2 (s2 - h), ai2 = -- 7T — r-r- ( 31 ) 

1 - {si,s 2 ) 

In this proof, we make a general position assumption that the spheres s and 
S12 do not touch and the spheres si and §2 do not touch. Write the linear 
dependence condition as 

§12 = As + fi§i + vs2- (32) 
Scalar product of this with s±, §2 leads to: 

1 + A = yu + v(si,S2) = — £t(si, §2) — v => n = -v = ——^— 

1 — \Sx,S2) 

Similarly, a scalar product of eq. (|3*2*|) with s, §12 leads to: 

— v = A — (s, §12) = 1 — X(s, §12) =>■ A = 1. 
This yields eq. (|3*T|) . □ 



6 Remarks on the continuous limit 

To perform a continuous limit in constructions of Discrete Differential Ge- 
ometry, one should think of the underlying lattice Z 2 of discrete surfaces 
as of (eZ) 2 and then send e — > 0. In such a limit it is common to assume 
that the discrete functions y(ni,ri2) on 1? approximate sufficiently smooth 
functions y{u\, u 2 ) on M. 2 , if rije = m. 

Moutard nets. It turns out, however, that the discrete Moutard equa- 
tion in the form with minus signs does not admit a nice continuous limit, 
while its close relative - discrete Moutard equation with plus signs - does: 

1/12 + y = 012(2/1 + V2)- (33) 

2 

Indeed, if a\ 2 = 1 + \q\ 2 , then the above equation is re-written as 

-2 (yi2 - vi - 2/2 + y) = -7012(2/12 + 2/1 + 2/2 + 2/), 
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which in the limit e — > tends to the Moutard equation |M] : 

did 2 y = qy. (34) 

Equations Q and (|33|) are actually related by a very simple transformation 
of dependent variables: y{n\,n 2 ) 1— > (— l) n2 y(ni, 712). As long as y play the 
role of homogeneous coordinates, such a transformation does not change the 
corresponding points in the projective space, therefore in the formulation of 
Theorem |21 we could use eq. (|33j) in place of eq. (JTJ), as well. However, 
eq. (|33jl cannot be posed on all faces of a three-dimensional lattice, and it 
is not multidimensionally consistent. As a consequence, the right way to 
construct discrete Moutard transformations for nets governed by eq. (|33|) is 
to consider multidimensional T-nets and then perform a change of variables 
y{n) >—¥ (— l) n ' 2 y(n) which breaks the symmetry among lattice directions. 
This leads to changing the signs in some of the formulas © for discrete 
Moutard transformations, which turn into 

yf - y = h(y + - yi), y£ + y = b 2 (y + + y 2 ). (35) 

An advantage of this break of symmetry is again the possibility to perform 
a continuous limit in eqs. (|35|) . recovering the formulas of the Moutard 
transformation for eq. (|34jl. see [Mj : 

d\y + + dxy = pi(y + - y), d 2 y + - d 2 y = p 2 (y + + y)- (36) 

Isothermic nets. A similar situation takes place for discrete isothermic 
nets. In order to enable the continuous limit to smooth isothermic surfaces, 
one should start with a discrete two-dimensional isothermic net with embed- 
ded elementary quadrilaterals. It is convenient to represent their negative 
cross-ratios as 

g(/,/i,/i2,/ 2 ) = -— , (37) 

a 2 

with positive labels a\ and a 2 . Formally, this means nothing more than 
changing the notation a 2 1— > — a 2 . But this operation puts some formulas 
into the form which allows for a continuous limit. If one keeps the both 
formulas 

l/i - /I 2 = aissi, I/2 - /| 2 = a 2 ss 2 , 

then this implies a slight modification in the definition of the function s, 
namely s(u) 1— > (— l) U2 s(u), which assures the positivity of s in case of 
positive ai, a 2 . Only upon this modification does the positive function 
s have a continuous limit. In this limit we recover a characterization of 
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isothermic surfaces as those curvature line parametrized surfaces for which 
the first fundamental form is conformal, possibly upon a re-parametrization 
of curvature lines: 



(dif,d 2 f) = 0, \d 1 f\ 2 = a 1 s 2 , \d 2 f\ 2 = a 2 s 2 , (38) 
with on = a.i{ui). Eq. (|2*2*|) for the dual surface turns into 

\hfV ssi \S 2 f\ 2 ss 2 

which is a direct discrete analogue of equations 

a f* dl f dl f a t* d2 f d2 f 

\dif\ 2 s 2 \d 2 f\ 2 s 2 

defining the dual isothermic surfaces in the smooth case. 

Remark. In the smooth case the functions «i, ot 2 in eq. (|38|) can be 
absorbed into a re-parametrization of the independent variables U{ 1— > y>i{ui) 
(i = 1,2), by which one can always achieve that ol\ = a 2 = 1, so that 
the first fundamental form of the surface / is conformal. Of course, such a 
re-parametrization is not possible in the discrete context. Nevertheless, one 
can consider a narrower class of discrete isothermic surfaces, characterized 
by eq. (|3*7j) with a± = a 2 = 1: 

?(/,/!, /l2,/ 2 ) = -l. (41) 

This condition (all elementary quadrilaterals of / are conformal squares) 
may be regarded as a discretization of the conformality of the first funda- 
mental form. Eq. (|41j). being a particular case of eq. (|2*Tj) with a special 
labelling, enjoys all the properties of the general case. However, it is im- 
portant to observe that it is not 3D consistent with itself, i.e., it cannot 
be imposed on all faces of a 3D cube. Indeed, if a\ja 2 = —1, then it is 
impossible to have additionally a 2 /a^ = —1 and a\/otz = —1. 

Note that the above change of signs yields also the similar modification 
in the lift s, which therefore satisfies the discrete Moutard equation with 
plus signs (|33[) . admitting the continuous limit. This characterization of 
isothermic surfaces is due to Darboux |Daj : 

Theorem 22 (Isothermic surfaces as Moutard nets in the light 
cone) A surface f : R 2 — > M. N is isothermic, if and only if there exists 
a function s : M 2 — > R such that the lift of f into the light cone 

S = a - 1 (/ + e + |/| 2 e oo ):R 2 -.L jV + 1 ' 1 
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satisfies a Moutard equation {3$ - Darboux transformations of an isothermic 
surface f are obtained by projection to 1^ from Moutard transformations 
of s within the light cone L +1,:L . 

L-isothermic surfaces. Analogously, we arrive at a (apparently new) 
characterization of smooth L-isothermic surfaces. 

Theorem 23 (L-isothermic surfaces as Moutard nets in the La- 
guerre quadric) A surface enveloping a two-parameter family of planes 
P : R 2 — > {planes in R 3 }, with P = {x G R 3 : (v,x) = d}, is L-isothermic, 
if and only if there exists a function p : R 2 — > R such that 

p = p -\ v + 2deoo + \v\e 6 ) : R 2 -» L 4 ' 2 

satisfies a Moutard equation \34\) - 
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